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We consider continuously monitored quantum systems and introduce definitions of work and heat 
along individual quantum trajectories that are valid for coherent superposition of energy eigenstates. 

We use these quantities to extend the first and second laws of stochastic thermodynamics to the 
quantum domain. We illustrate our results with the case of a weakly measured driven two-level 
system and show how to distinguish between quantum work and heat contributions. We finally 
employ quantum feedback control to suppress detector backaction and determine the work statistics. 
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Thermodynamics is, at its heart, a theory of work and 
heat. The first law is based on the realization that both 
quantities are two forms of energy and that their sum 
is conserved. At the same time, the fact that entropy, 
defined as the ratio of reversible heat and temperature, 
can only increase in an isolated system is an expression 
of the second law mm- In classical thermodynamics, 
work is defined as the change of internal energy in an 
isolated system, W = At/, while heat is introduced as 
the difference, Q = At/ — W, in a nonisolated system. 
Thermal isolation is thus crucial to distinguish W from 
Q. In the last decades, stochastic thermodynamics has 
successfully extended the definitions of work and heat 
to the level of single trajectories of microscopic systems 
[^. In this regime, thermal fluctuations are no longer 
negligible and the laws of thermodynamics have to be 
adapted to fully include them. The second law has, for 
instance, been generalized in the form of fluctuation the¬ 
orems that quantify the occurrence of negative entropy 
production [3]. A particular example is the Jarzynski 
equality, (exp(—/3IT)) = exp(—/3AF), that allows the 
determination of equilibrium free energy differences AF 
from the nonequilibrium work statistics in systems at ini¬ 
tial inverse temperature (3 [S]. The laws of stochastic 
thermodynamics have been verified in a large number of 
different experiments, see Refs. m H] and the review O. 

The current challenge is to extend the principles of 
thermodynamics to include quantum effects which are 
expected to dominate at smaller scales and colder tem¬ 
peratures. Some of the unsolved key issues concern the 
correct definition of quantum work and heat, means to 
distinguish between the two quantities owing to the blur¬ 
ring effect of quantum fluctuations, and the proper clar¬ 
ification of the role of quantum coherence. A variety of 
approaches have been suggested to tackle these problems 
[M], and quantum work statistics has been measured 
in isolated systems in two pioneering experiments using 
NMR and trapped ions [52]. A new approach may 
emerge from the possibility of weakly monitoring quan¬ 
tum systems. Recently, individual quantum trajectories 
of a superconducting qubit in a microwave cavity have 


been observed using weak measurements |23l |24| . These 
measurements only slightly disturb quantum systems ow¬ 
ing to the weak coupling to the measuring device [25] . 
They hence allow to gain information about states with¬ 
out projecting them into eigenstates. They have been 
successfully employed to explain quantum paradoxes [26] , 
detect and amplify weak signals [IZI [55], determine a 
quantum virtual state |29j . as well as directly measure 
a wave function |30j . Motivated by the two experiments 
[551151] , we here investigate the first and second law for 
continuously monitored quantum systems and aim at de¬ 
veloping a quantum stochastic thermodynamics based on 
quantum trajectories. Such an extension faces several 
technical difficulties. First, since a weakly measured sys¬ 
tem can be in a coherent superposition of energy eigen¬ 
states, energy is not a well-defined concept along a single 
quantum trajectory. Furthermore, even in the absence of 
an external environment, a continuously monitored quan¬ 
tum system is not isolated and the detector backaction, 
albeit small, will perturb its dynamics [5^. As a result, 
its time evolution will be nonunitary and energy, in the 
form of heat, will be exchanged with the detector. 

In the following, we introduce suitable and consistent 
definitions of work and heat contributions to the quan¬ 
tum stochastic evolution of a weakly measured system 
that is externally driven. We use these definitions to de¬ 
termine the distributions of quantum work and heat for 
a two-level system, and demonstrate the general validity 
of the Jarzynski equality, hence of the second law. We 
finally use the tools of quantum feedback control [3T] to 
suppress detector backaction and thus effectively achieve 
thermal isolation of the system. This provides a practi¬ 
cal scheme to experimentally test our definitions of work 
and heat along individual quantum trajectories. 

Quantum work and heat. We consider a system with 
time-dependent Hamiltonian Fit that is initially in a ther¬ 
mal state at inverse temperature /3, po = exp(— 
where Zq is the partition function. The system is driven 
by an external parameter At during a time r. At the en¬ 
semble level, quantum work and heat are introduced by 
considering an infinitesimal variation of the mean energy. 
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U = (H) = tvlptHt} ^133]: 

dU = tr {ptdHt} + tr {dptHt} = 5W + 5Q. (1) 

Heat is further related to entropy S = —kir{pt Inpt} via 
5Q = TdS |32l |33j. For an isolated system with uni¬ 
tary dynamics heat vanishes, since dS = 0, and therefore 
dU = 5W in agreement with classical thermodynamics 
mm- Heat therefore appears to be fundamentally asso¬ 
ciated with the nonunitary part of the dynamics. 

At the level of individual realizations, energy is a 
stochastic quantity owing to thermal and quantum fluc¬ 
tuations. The distribution p(u) of the total energy change 
u may be determined by performing projective measure¬ 
ments n„ and Hm, with outcomes E° and at the 
beginning and at the end of the driving protocol [31 [M] , 

p{u) = Y, PknPl5{u - AF;^.„). (2) 

m,n 

Here = tr {n„po} denotes the probability of the eigen¬ 
value E^, „ = trjHmPn.T} the transition probability 

from state n to m, with pn^r the time evolved projected 
density operator p^p = n„pon„/P°, and AEm,n = 
E^ — E^ the energy difference. For unitary dynamics, 
p(u) reduces to the work distribution p{W), but, in gen¬ 
eral, Eq. ([^ does not allow to distinguish work from heat. 
In the following, we generalize Eq. ([^ and identify work 
and heat for a weakly measured system. 

A quantum system continuously monitored by a quan¬ 
tum limited detector may be assigned, for each individual 
trajectory, a conditional density operator pt that reduces 
to the usual density operator pt when averaged over all 
the trajectories, pt = {{pt}) [311135]. The evolution of pt 
is commonly described by a stochastic master equation 
that contains a random parameter ^ (t) that accounts for 
the detector shot noise, see Eqs. (l9|-@ below for an 
example. An important observation is that such master 
equation has a unitary component, corresponding to the 
dynamics generated by the system’s Hamiltonian, and a 
nonunitary part that stems from the continuous coupling 
to the detector. For an infinitesimal time step, these two 
contributions are additive and may be written as, 

dpt = SW[pt]dt + 6Q[pt]dt, (3) 

where (5W[;Ot] and (5Q[pt] are operators associated with 
the respective unitary and nonunitary parts of the dy¬ 
namics. We identify them as corresponding to work and 
heat at the level of an infinitesimal quantum trajectory. 
This separation cannot be directly extended to the entire 
(time integrated) trajectory, since the stochastic master 
equation is generally a nonlinear function of the operator 
Pt- However, when averaged over quantum fluctuations, 
Eq. § allows to extend the first law Q to single real¬ 


izations of the stochastic measurement outcome, 

dilt =tr {Ht{pt-dt + dpt)} - tr {Ht-dtPt-dt} 

=tr {pt-dtdHt} + tr {HtSWdt} + tr {HtS^dt} 

=5Wt+5Qu (4) 

where in the second line dPdt = Ht — Ht-dt and the mid¬ 
dle term tr{Ht5Wdt} = 0 since SW is unitary [36|. In the 
last line, dWt = ti {pt-dtdHt} and SQt = tr{E[tSQdt}, 
indicating that work is related to a change of the Hamil¬ 
tonian, as expected, and heat to the nonunitary <5Q. 
Equation Q is a direct extension of stochastic thermo¬ 
dynamics to the quantum domain. The first law Q is 
recovered when Eq. (§ is averaged over both stochastic 
and quantum fluctuations. The integrated work and heat 
contributions to the changes in transition probabilities, 
dPm,,n = Pm,n “ Pm,n^ ^ith „ = 6n,m the initial tran¬ 
sition probability, can be further obtained from Eq. ([^ 
by carefully adding all the different terms (see Ref. [37[T. 
We find, for each individual quantum trajectory, 

dP^,r^ = 6PZn+SP2,u, (5) 

with the two quantities. 
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f dt<5W[pt]| 
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= tr <j 


f (iti5Q[pt]| 


These expressions depend explicitly on the quantum tra¬ 
jectory Pt which we stress by using the notation 6 instead 
of d. They provide an unambiguous way to distinguish 
between work and heat at the level of a single trajectory. 
Remarkably, they are valid even if the system remains in 
a coherent superposition of energy eigenstates, that is, 
when its energy is ill-defined. Equation ([^ holds for the 
trajectory averaged quantities dPm,n = dP^^ + dP^„, 
with dP^ ,^ = ((^Rm_„)), a = W,Q. This averaged dis¬ 
tinction between work and heat contributions to tran¬ 
sition probabilities requires access to single trajectories, 
thus cannot be established directly at the ensemble level. 

The second law in the form of the Jarzynski equal¬ 
ity, f dlFp(IF) exp(—/3VF) = exp(—/3AF), immediately 
follows from Eq. (|^ for an isolated system [Sj. How¬ 
ever, the equality is not satisfied for an open system 
with nonunitary dynamics owing to the heat term [34] . 
The second law may be restored by replacing P^ „ by 
PZn = P^,n + ^PZn^ by Setting „ to zero 

at each time step, see Fig. 3. We next show how quantum 
work and heat may be identified theoretically, by numer¬ 
ically analyzing a weakly measured two-level system, and 
experimentally, by means of quantum feedback control. 

Application to a monitored qubit. In order to illus¬ 
trate our approach, we consider a driven two-level sys¬ 
tem S with Hamiltonian, Ht = euz + \t<^x-, where At is 
the external driving and at the usual Pauli matrices. The 
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FIG. 1. (color online) First law for a weakly measured qubit, 
a) Infinitesimal change of work, heat and energy along a single 
quantum trajectory pt \ for each realization dUt = SWt -\-5Qt, 
Eq. (4). b) Corresponding signal/(f) in the detector. Param¬ 
eters are S'o/A/^ = 2.5 • 10®df, h/g = 1.6 • lO^df, Yijt = 10^dt, 
u = 8 and r = 3 • 10^dt (see main text). 


system is continuously coupled to a quantum limited de¬ 
tector D via the interaction Hamiltonian Hf. 

H = Ht + HD + a,Hi, ( 8 ) 

where, without loss of generality, we identify as the 
system’s observable that is monitored by the detector. 
The effect of the detector is fully characterized by the 
averaged signals, (Ii, I 2 ), and Gaussian noises, {Si, S 2 ), 
measured when the qubit is in the two eigenstates, (|1), 

12)), of the measured observable. We assume to be in 
the weak measurement regime, i.e. at time scales smaller 
than the measurement time tm = {Si + S 2 )/{Ii — ^ 2 )^. 
For concreteness and simplicity, we will interpret the 
qubit in Hamiltonian ([^ as describing a double quan¬ 
tum dot sharing a single electron and interacting with 
a quantum point contact (QPC), but it can also be ap¬ 
plied to a qubit coupled to a microwave resonator [55] in 
a circuit QED set-up as in the experiments (251121]. We 
accordingly identify the configurations where the elec¬ 
tron occupies only one dot by {az) = ±1. Coherent 
superpositions of the two are possible. The detector 
monitoring the occupation of the dots is a voltage bi¬ 
ased QPC with Hamiltonian [55H55] Ho = Eiajai + 
Eralar + J2i r +«; Or) and the interaction term 

reads Hi = ^5Vt/2{a\.ai ajar). The signal in the 

detector is the current /(f) across the QPC, with aver¬ 
ages /i( 2 ) = 2TTn\^^)Pipre^V/h = e^ri( 2 )H//i and noises 
5 'i( 2) = e(l — '7i(2))di(2). Here ppr are the densities of 
states in the left and right electrodes, and Ti( 2 ) are the 
dimensionless transmission probabilities across the QPC. 

Under the assumption of a weakly coupled detector, 
the detector signal is a random variable, and the evo¬ 


lution of the system depends on the specific realization 
of the stochastic process. This is captured by a well- 
established Bayesian formalism [HI [42] which describes 
the evolution of the system conditional to the detector’s 
outcome in terms of a nonlinear stochastic differential 
equation for the system’s density matrix p{t). In the Ito 
formulation, we have [iiliig, 


hi = - 2^ Im {P12) + Pii(l - ( 9 ) 

h So 

. .A(f) , (A/)2 

P 12 —2i-pi2 — z—^ (I — 2pii) — Pi2—ro — 
h h 4c»o 

+ {I - 2hi)h2-^^{t), ( 10 ) 

Oq 


where A/ = I 2 — h and ^(f) is the white noise of the 
detector’s signal with (^(f)) = 0, = cr^(5(f — f') 

and (T = ^ySo/2. The detector current /(f) is further, 

= do + —(2pii — I)-I-^(f). (II) 

For each realization of the measurement outcome, 
Eqs. ([^ and (101 allow to identify the unitary and 
nonunitary contributions to the time evolution, since 
the nonunitary part is proportional to A/. We rewrite 
Eq. (|^ as dh = (5W[/5t](ff -I- (A//5'o)5M[pt](if, and iden¬ 
tify SW with the work done by the driving At along an 
infinitesimal trajectory and (5Q = (A//S'o)(5M as the 
heat associated with the detector backaction of the de¬ 
tector. Due to the nonlinearity of the stochastic mas¬ 
ter equation, we can only determine the distributions of 
work and heat numerically. We specify the driving as 
At = g(l/cosh(jz(I — f/r)), where r is the duration of 
the experiment, and reformulate equations and ( [To] ) 
in the Stratonovich form im \M- We solve them nu¬ 
merically by the Monte-Carlo method for an ensemble of 
300 realizations of the random signal /(f) in the interval 
f/r S [0,1] using a time step dt = O.OI. The results for 
work, heat and energy along a given quantum trajectory, 
Eq. ([^, are shown in Fig.[^ while those for the work and 
heat contributions to the transition amplitudes, Eqs. (j^- 
([ 7 ]), are presented in Fig. (see Ref. [57] for details). 

Figure [^) demonstrates the reconstruction of quan¬ 
tum averaged work and heat changes, SWt and 
along a single quantum trajectory, based on the defini¬ 
tions given in Eq. (4). The corresponding signal /(f) 
in the detector is displayed in Fig. &)■ Contrary to 
the case of an isolated system for which dll = 5W, the 
heat contribution SQt is here clearly visible. Equation([^ 
holds for each individual realization and thus extends the 
first law of stochastic thermodynamics to the quantum 
regime. Figure]^) shows the unambiguous distinction of 
the work and heat contributions, SP^.^ and SP^ .^, eval¬ 
uated via Eqs. to the final transition probability 

Pm n- We stress that, although P{^ „ is always positive, 
as a proper probability should be, the work and heat con¬ 
tributions need not be: the probability to go from state 
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FIG. 2. (color online) a) Averaged final transition probabil¬ 
ities (yellow) for a continuously monitored qubit with 

their work and heat contributions, (blue) and „ 

(red), and the initial transition probability, P^^„ = Smn (pur¬ 
ple). The first law like equation dPm,n = SP^„ + SP^ .^ is 
verified, b) Work and heat contributions, 5P^^^ and SP^ „ 
at the single trajectory level. Same parameters as in Fig. [T] 



FIG. 3. (color online) Transition probabilities P^n for the 
weakly measured qubit with (orange) and without (yelllow) 
feedback control for a) ri = 1.4 • lO^dt and b) r 2 = 2.5 • lO^dt. 
The feedback strength is / = 3. The isolated, unitary, case 
(red) is shown as a reference. The feedback loop effectively 
suppresses the detector backaction and the associated heat 
exchange, achieving thermal isolation. 


n to m at time r can, for instance, be smaller than the 
initial transition probability |43j . Note that a quantity, 
dPZ,n = P^Ja - Pm%^OL = W,Q, that only depends on 
initial and final times, cannot be defined, reflecting the 
fact that there are no heat or work operators. 

Quantum feedback control. In classical thermody¬ 
namics work is associated with the variation of the in¬ 
ternal energy of the isolated system mm- After having 
shown above how heat can be theoretically identified, we 
next take advantage of a feedback loop protocol to sup¬ 
press the detector backaction [3T], offering a scheme to 
reach isolation experimentally. Quantum feedback has 
recently been demonstrated experimentally for a super¬ 


conducting qubit [IS]. Specifically we control the ampli¬ 
tude, g, of the system’s driving depending on the contin¬ 
ues detector outcome, i.e. g ^ gt = {1 — fA(pt)g, where 
/ is the feedback strength and Aipt the phase difference 
between the actual vector (with backaction) and desired 
vector (without backaction) in the Bloch sphere of the 
qubit (see Refs. Eznn] for details). This allows to op¬ 
erationally counter the effects induced by the continuous 
monitoring. From a thermodynamic point of view, the 
feedback adds an extra amount of work that exactly can¬ 
cels the heat contribution to the transition probabilities. 

Figure shows the numerically simulated final tran¬ 
sition probabilities „ for the weakly measured qubit 
with (orange) and without (yellow) quantum feedback 
for two driving times. We observe in both cases that the 
feedback process effectively suppresses the heat contribu¬ 
tions (identified in Fig. and that the transition proba¬ 
bilities agree with those of the isolated system with uni¬ 
tary dynamics (brown). Quantum feedback control thus 
appears as a powerful tool to determine the statistics of 
the work done by the external driving in a continuously 
monitored system. The heat statistics can be further eas¬ 
ily obtained by measuring the undriven system, that is, 
when no work is performed and dUt = SQt. 

The above findings can be directly used to verify the 
quantum Jarzynski equality for the driven qubit. Since 
any measurement induced heating is prevented by the 
feedback, only the initial inverse temperature /3 of the 
system matters. Determining the quantum work statis¬ 
tics via Eq. (2), we find AFi = —0.488 and AF 2 = 
—0.496 with feedback control for ti = 1.4 • lO^dt and 
T 2 = 2.5 • lO^dt and AF = —0.495 in the unitary case (for 
/3 = 10). The excellent agreement demonstrates the cor¬ 
rectness of the definitions of work and heat, and confirms 
the second law for a weakly measured quantum system. 

Conclusions. We have extended the laws of stochas¬ 
tic thermodynamics along individual quantum trajecto¬ 
ries of a weakly measured system. We have shown how 
to distinguish work and heat contributions to both the 
energy changes and the transition probabilities. We have 
further demonstrated the usefulness of our approach with 
the analysis of a driven qubit and introduced methods to 
identify work from heat numerically as well as experi¬ 
mentally with the help of quantum feedback control. 
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Note added. While completing this manuscript, we 
became aware of a preprint by Elouard et al. [46] that 
also discusses quantum stochastic thermodynamics. 
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